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Abstract 



The Newtonian Euler-Poisson equations with attractive forces are the classical models for 
the evolution of gaseous stars and galaxies in astrophysics. In this paper, we use the integration 
method to study the blowup problem of the iV-dimensional system with adiabatic exponent 
7 > 1, in radial symmetry. We could show that the C 1 non-trivial classical solutions (p, V), 
with compact support in [0, i2] , where R > is a positive constant with p(t,r) = and 
V(t, r) = for r > R, under the initial condition 



with an arbitrary constant n > max(iV — 2, 0), 

blow up before a finite time T for pressureless fluids or 7 > 1. Our results could fill some 
gaps about the blowup phenomena to the classical C solutions of that attractive system with 
pressure under the first boundary condition. 

In addition, the corresponding result for the repulsive systems is also provided. Here our 

result fully covers the previous case for n = 1 in "M.W. Yuen, Blowup for the Euler and 
*E-mail address: nevctsyuen@hotmail.com 
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1 Introduction 

The compressible isentropic Euler (6 = 0) or Euler-Poisson (S = ±1) equations can be written in 
the following form: 

p t +V ■ (pu) =0 

S p[u t + (u ■ V)u]+VP =pV$ ( 2 ) 
A$(t,x) = 5a{N)p 

where a(N) is a constant related to the unit ball in R N : : a(l) = 1; cv(2) = 2-7T and for TV > 3, 

a(N) = N(N - 2)Vol{N) = N(N - 2) + , (3) 

where Vol(N) is the volume of the unit ball in R N and T is a Gamma function.. As usual, 
p = p(t,x) > and u = u(t,x) £ H, N are the density and the velocity respectively. P = P(p) is 
the pressure function. The 7-law for the pressure term P{p) could be applied: 

P(p)=Kp^ (4) 

which the constant 7 > 1. If K > 0, we call the system with pressure; if K = 0, we call it 
pressureless. 

When S = — 1, the system is self-attractive. The equations @ are the Newtonian descriptions of 
gaseous stars or a galaxy in astrophysics [5] and [3]. When 6 = 1, the system is the compressible 
Euler-Poisson equations with repulsive forces. It can be used as a semiconductor model [B]. For 
the compressible Euler equation with 6 = 0, it is a standard model in fluid mechanics |16) . And 
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the Poisson equation @3 could be solved by 



$(t,x)=6 G(x-y)p(t,y)dy (5) 

JR N 



where G is Green's function: 

xl N = l 



G(x) = { 



Here, the solutions in radial symmetry could be: 



log |4 N = 2 (6) 



with the radius r = QZiLi x 1 

The Poisson equation $Z§ 3 becomes 



p = p(t,r) and u = -V(t,r) =: -V (7) 
r r 

1/2 



r^ 1 ^ (t, a:) + (N - 1) r^-^^a (JV) Jpr^ 1 (8) 
a (iV) <5 



$ r = 



r JV- 



^ [ r p(t,s)s N - 1 ds. (9) 
Jo 



By standard computation, the systems in radial symmetry can be rewritten in the following form: 

N - 1 

Pt + Vp r + pV r H pV = 

r (10) 

+ w r ) + p r Go) = p$ r (p). 

In literature for constructing analytical solutions for these systems, interested readers could 
refer to [T3] , Q2] , [23] ; [H] and [M] • The local existence for the systems can be found in [TB] , [T8] , 
Q] and [12]. The analysis of stabilities for the systems may be referred to [11], [19], [20], [21], [9], 
dm, [23], 0, H and |25] . 

In literature for showing blowup results for the solutions of these systems, Makino, Ukai and 
Kawashima firstly defined the tame solutions [19 for outside the compact of the solutions 

V t + VV r = 0. (11) 

After this, Makino and Perthame continued the blowup studies of the "tame" solutions for the Euler 
system with gravitational forces |20j . Then Perthame proved the blowup results for 3-dimensional 
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pressureless system with repulsive forces [2T] (5=1). In fact, all these results rely on the solutions 
with radial symmetry: 

V t + VV r =^^ f p(t,s)s N - 1 ds. (12) 
r Jo 

And the Emden ordinary differential equations were deduced on the boundary point of the solutions 
with compact support: 

= J^T< fi (°< R o) = Ro>0, R(0, R ) = (13) 

where 4? := V and M is the mass of the solutions, along the characteristic curve. They showed 
the blowup results for the C 1 solutions of the system (|10l) . 

In 2008 and 2009, Chae, Tadmor and Cheng in [3| and [7] showed the finite time blowup, for the 
pressureless Euler-Poisson equations with attractive forces (S = —1), under the initial condition, 

S :={ae R N \ p (a) > 0, O (a) = 0, V-it(0,x(0) < 0} ^ (14) 



where fl is the rescaled vorticity matrix (£l ij) = ^(9;Uq— 9jUg) with the notation u = (u 1 , u 2 , u N 
in their paper and some point a:o- 

They used the analysis of spectral dynamics to show the Riccati differential inequality, 

Ddivu 1 2 

< (dxvu) . (15) 

Dt ~ N y ' K ' 

The C 2 solution for the inequality dTSl blows up on or before T = —NJ (V • «(0, xo(0)). But, their 
method cannot be applied to the system with the system with pressure or the classical C 1 solutions 
for the system without pressure. 

We are in particular interested in the instabilities of the Euler-Poisson equations in R 3 . For the 
instabilities of the global classical solutions for 7 > 4/3 is shown by the second inertia functional 

H{t) = I p(t,x) \x\ 2 dx (16) 

with the constant total energy 

E = Eit) = J Qp M 2 + + ip$) dx (17) 
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[20] and [5]. For the critical case 7 = 4/3, any small perturbation which makes the energy positive 
would cause that the solutions go to positive infinity, which implies an instability of such a state. 
And the instability for the stationary solution for 7 = 6/5 is obtained |14j . But the nonlinear 
stability result for the other cases is unknown. It was shown that any stationary solution is stable 
for 7 € (4/3,2) and unstable for 7 G (1,4/3) in linearized case [5] and |15j . 

On the other hand, Yuen [55] in 2010 used the integration method to show that with the initial 
velocity 



the solutions with compact support to the Euler (5 = 0) or Euler-Poisson equations with repulsive 
forces (6 = 1) blow up before the finite time T. 

In this article, we could further apply the integration method in [25] for the system with 
attractive forces (S = —1) to fill some gaps about the blowup phenomena for the classical C 1 
solutions of the system with the pressure 7 > 1 or pressureless fluids: 

Theorem 1 Consider the Euler-Poisson equations with attractive forces (S = —1) S10\) in R N . 
The non-trivial C 1 classical solutions (p, V), in radial symmetry, with compact support in [0,R], 
where R > is a positive constant (p(t,r) — and V(t,r) — for r > R) and the initial velocity: 




(18) 




(19) 



with an arbitrary constant n > max(N — 2, 0) and M is the total mass of the fluids, 



blow up before a finite time T for pressureless fluids (K = 0) or 7 > 1. 



Here, the boundary condition for the fluids: 



p(t, r) = and V(t, r) = for r > R 



(20) 



is called non-slip boundary condition [8] and [5]. 
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2 Integration Method 



We just modify the integration method which was initially designed for the Euler-Poisson system 
with repulsive force (S = 1) in [25] to obtain the corresponding blow results for the ones with 
attractive forces (S = —1). 

Proof. We only use the density function p(t,x(t;x)) to preserve its non-negative nature as we 
take integration for the mass equation ([2]) i : 

Dp 



Dt 



pV • u = 



(21) 



with the material derivative, 



D__d_ 

Dt~ dt 



+ (u • V) 



to have: 



p(t, x) = po(x o (0, x )) exp ( - J V • u(t, x(t; 0, x ))dt ) > 



(22) 



(23) 



for pq(xq(0, Xq)) > 0, along the characteristic curve. 

Then we could control the momentum equation (|10[1? for the non-trivial solutions in radial 
symmetry, po ^ 0, to obtain: 

V t + VV r + K 1 p~ 1 ~ 2 Pr = $ r 
d ,1. 



v t + ^qv 2 )+K 7 p^- 2 p r = ^ r 

r n V t + r n ^-(\v 2 ) + K lr n p^- 2 p r = r"$ r 
or 2 



(24) 
(25) 
(26) 



with multiplying the function r" with n > 0, on the both sides. 

We could take the integration with respect to r, to equation (|2l)|) . for 7 > 1 or K > 0: 

/ r n V t dr+ r n — (-V 2 )+ K lr n p^ 2 p r dr = r n t> r dr 
Jo Jo dr 2 Jo Jo 



r n V f dr- 



d ,1 



dr 2 



Gv 2 ) 



7-1 



dp 1 



r r 



a{N)r n [ r 



r N-l 



p{t 1 s)s N - 1 ds 



(27) 
dr (28) 



with the estimation for the right hand side of ([28]) : 



r r 



a(N)r n f r 



r N — l 



p(t,s)s"- l ds 



dr < 



a{N)r n 

r N-l 



p(t,s)s N - 1 ds 



dr = I r 
'0 



R nn-N+2i\/r 

"- N +Hldr R M - 



n-N + 2 



(29) 
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where M is the total mass of the fluids and the constant n > max(7V — 2, 0), 



a(N)r r ' 

~.N-1 



p(t, s)s"- i ds 



dr > 



Rn - N+ 2 M 
' n-N + 2 



to have 



( R r n V t dr+ f 
Jo Jo 



r " v "' r+ 1 rn i { \ v2) 



R Kjr n „ , R n ~ N+2 M 
-dp 1 ' 1 > -- 



(30) 



(31) 



Jo 7-1"" " n-N + 2 
We notice that this is the critical step in this paper to obtain the corresponding blowup results for 
the Euler system with attractive forces. 

Then, the below equation could be showed by integration by part: 

j R r n V t dr - ^J Q R nr n - 1 V 2 dr + \ [R n V 2 (t,R) - 0™ • V 2 (t,0)] 
Io R ^r^P^dr + ^ [Trpi-\t, R) 0" ■ pr-^t, 0)] > -^+2". 
The above inequality with the boundary condition with a uniform compact support (V(t, R) = 
and p(t,R) = 0), becomes 



(32) 



M 



[ r n 
Jo 



V t dr 



1 



R 



nr n - x V 2 dr - 



R 



K-ynr 71 ' 1 ^ , , R n - N+2 M 
p < dr > 



7-1 



n-N + 2 



(33) 



d 1 



d 
~dt 

,-R 



— I r n Vdr - - I nr n ~ x V 2 dr - 



K^nr 11 - 1 „ , , R n - N+2 M , , 

-pi- x dr > — - (34) 



- / Vdr n+1 - - [ - 11 - V 2 dr n+1 + 
lio 2i (n + l)r 



7-1 

R n-N+2 M 



n-N + 2 



dtn + 

for n > and 7 > 1 or K = 0. 

For non-trivial initial density functions p > 0, we have: 



n-N + 2 



> 



f 

Jo 



Kjnr 



n-l 



7-1 



-p r '- 1 dr>0 (35) 



d 1 
dt n + 1 J 



Vdr 



n+l 



1 



2 Jo (n+l)r 



nn-N+2iif 



n-N + 2 



d 1 
dtn+l Jo 



Vdr 



R n-N+2 M ,fl 

J > / 



n- TV + 2 ~ 7 2(n+l)r 



R n-N+2 M 



V 2 dr n+1 > 



2{n+l)RJo 



-R 



V 2 dr r 



Vdr 



n+l 



n-N + 2 ~ 2R 



f 

Jo 



We could denote 



£T := H(t) = [ r n Vdr = — !— [ Vdr n 
Jo n+l J 

and apply the Cauchy-Schwarz inequality: 



(36) 

(37) 
(38) 

(39) 



/' 

Jo 



V ■ ldr n+1 



< 



j\ 2 dr n+ ^j (^J\dr n+1 ^ 



1/2 



(40) 
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ldr n+1 




J R V 2 dr n+1 ^j 









1/2 



R 

I I V 2 dr n+1 



1/2 



-IN V2 



R n+1 ~ Jo 

n{n + l) 2 H 2 <JL f R v2drn+1 



2R n + 2 ~ 2R J a 



for letting equation (f38|) to be 



R 



d , . R n - N+2 M n . 



r2 ^„ +1 ^ n(n + l) 2 g 2 



with inequality (Hi)) . 



d n(n + l)H 2 R n ~ N+2 M 
dt ~ 2R n + 2 n-N + 2' 



If we require the initial condition 



H ° > V n(n + l)(n-7V + 2) : 



(41) 

(42) 

(43) 
(44) 



(45) 



(46) 



(47) 



it is well-known for that the solutions for the Riccati differential inequality (|4"6l blow up before a 
finite time T: 

Km H(t) = +oo. (48) 

t->T- 

Therefore, we could show that the C 1 solutions blow up before a finite time T. 
This completes the proof. ■ 

Additionally, it is clear to see that we could further apply the integration method to extend 
Yuen's result [5] as the following theorem: 



Theorem 2 Consider the Euler (8 = 0) or Euler-F 'oisson equations with repulsive forces (8 = 1) 
(0) in R . The non-trivial classical solutions (p, V) , in radial symmetry, with compact support in 
[0,-R], where R > is a positive constant (p(t,r) = and V(t,r) = for r > R) and the initial 
velocity: 



H = f r n V 
Jo 



dr>0 



(49) 
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with an arbitrary constant n > 0, 

blow up on or before the finite time T — 2R n+2 / (n(n + 1)Hq), for pressureless fluids (K = 0) or 
7 > 1. 

Proof. We just reverse the estimation for the equation (|2"8|) with repulsive forces (S = 1) to have 



p(t,s)s N -'ds 



dr (50) 



f r d 1 f K^r n 

\ r n V t dr+ r ' l -(-V 2 )+ — 7 — dp^ 1 > (51) 
Jo Jo dr 2 Jo 7-1 

for 5 > 0. 

For non-trivial initial density functions po > 0, we may obtain the corresponding result: 

dt ~ 2R n + 2 y ' 

> -2R n + 2 H 

~ n(n + l)H t - 2R n + 2 ' 1 ' 

Then, we could require the initial condition 

H = / r n V dr > (54) 
Jo 

for showing that the solutions blow up on or before the finite time T — 2R n+2 / (n(n + 1)Hq). 
This completes the proof. ■ 

We notice that Theorem 2 in this paper fully cover the previous case for n = 1 in [25] . Further 
researches are needed to have the corresponding results for the non-radial symmetric cases. 

Remark 3 // the global C 1 solutions with compact support whose radii expand unboundedly when 
time goes to infinity, the results in this paper could not offer the information about this case. 

Remark 4 The blowup results in this paper imply that the classical Euler-Poisson equations even 
for the ones with attractive forces could not be used to be a good modelling for the evolutions of 
gaseous stars with the relativistically large functional Hq which compares with the total mass of 
the fluid. Alternatively, the relativistic Euler-Poisson equations may be adopted for these cases to 
prevent these blowup phenomena. 
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On the other hand, the author conjectures that there exists a variational version for showing 
blowup phenomena in non-radially symmetrical cases with the first boundary condition. But, an- 
other nice functional form would be required to be designed to handle the corresponding problems. 
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